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1 The framework 



Kinetic equations are a widely used and highly successful method to understand 
dynamical aspects of quantum many particle systems. As a standard reference I 
mention the book by Ziman [1] which in large parts is based on the kinetic method. 

My article deals with the issue to derive such an equation in a well-defined limit 
from quantum dynamics. For sake of a simple presentation I restrict myself to a 
weak interaction potential. There is a parallel discussion for low density for which 
I have to point to the literature [21 13] . Let me first describe the general framework 
in admittedly rather vague terms. The broad picture will be made precise in the 
remainder of the article, where also the border between proven and conjectured will 
be traced out. 

We consider (nonrelativistic) quantum particles moving in three dimensions and 
interacting through a short range potential. The hamiltonian is of the generic form 



where H$ denotes the free hamiltonian. H\ is the interaction which is multiplied 
by a dimensionless strength parameter yfe. We want to study a system of many 
particles and, first, consider a spatially homogeneous situation, which is idealized 
as an infinitely extended system in a spatially translation invariant state. Amongst 
the many conceivable states there is a special class of states characterized by being 
quasifree, translation invariant, and time-invariant under the dynamics generated 
by H . Such a state is uniquely characterized by a mass density p, p > 0, and a 
momentum distribution W(k)/p, J dkW(k) = p. 

Under the full dynamics generated by H such a quasifree state changes in time 
but remains translation invariant. The change is slow, since the interaction is weak. 
By definition, the kinetic time scale is the shortest scale over which the change in the 
momentum distribution is of order 1. Since transition rates are roughly proportional 
to \^/~eHi\ 2 , the kinetic time scale should be e~ 1 t, in microscopic time units. On the 
kinetic time scale the momentum distribution evolves according to an autonomous 
nonlinear transport equation. Generically, on the basis of this equation, the long 
time limit will be the thermal equilibrium momentum distribution for Hq, i.e. for 
the noninteracting particles. Moreover, the kinetic equation makes a quantitative 
prediction on how thermal equilibrium is approached. 

At this stage the difficulty of a proof becomes already apparent. Since £ < 1, 
Equation ( 11. ip looks like a standard perturbation problem. But this perturbation 
has to be controlled over the very long times e~ 1 t. 

What happens if the initial state is not quasifree but still translation invariant? 
This situation is analysed in [I] for a weakly interacting Fermi gas on a lattice. If 
the initial state has good spatial clustering, then, under the i7 -dynamics for long 
microscopic times, the state becomes quasifree without yet changing the momentum 
distribution. The resulting dynamical picture is rather compelling: For "reasonable" 
initial states, the time-evolved state becomes rapidly quasifree. At kinetic times the 
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interaction makes itself felt, changing the momentum distribution but maintaining 
the quasifree property. At even longer times the system presumably tries to establish 
the equilibrium state for H, which is not strictly quasifree. For fixed small e, phys- 
ically the time scales are not perfectly separated. But mathematically a separation 
can be achieved through taking limits. The first regime is e = and t — > oo. The 
kinetic regime is e — > and t — > oo, et = t fixed, while the third regime is e fixed 
and small, while t — > oo. Needless to remark that the third regime is inaccessible 
mathematically. 

Usually, one also wants to include spatial variations, which however are slow 
on the scale of the typical distance between particles. For such slow variations, the 
i/o-dynamics can be approximated by semiclassics, which means that particles move 
independently with classical kinetic energy tv(k), u being the dispersion relation for 
the quantum particles. The kinetic time scale has been determined already. To have 
the if -dynamics of the same order as the collisions requires then a spatial variation 
of order e™ 1 in microscopic units. The momentum distribution varies now spatially 
on the kinetic scale. Thus W(k) is to be replaced by W(r,k). The collision term 
is strictly local in r, while the kinetic equation picks up the flow term — Vfcw(A;) • 
V r W(r, k) from the semiclassical approximation. Close to r the microscopic state 
is still quasifree and translation invariant, at least in approximation. 

2 Model hamiltonians 

There is some freedom in the choice of the microscopic model. Firstly, the position 
space is either M. 3 or some Bravais lattice which we take to be Z 3 . (We will comment 
on the spatial dimension below.) The average number of particles is finite but 
diverges as e~ 3 in the kinetic limit. Thus the Hilbert space for the system is the 
boson Fock space JF + , resp. fermion Fock space JF_, over either L 2 (M 3 ) or Z 2 (Z 3 ). 
Let a(x), a(x)*, x G M 3 , resp. x G Z 3 , denote annihilation and creation operators for 
the Bose, resp. Fermi, fields. Setting h — 1, they satisfy the commutation relations 

[a(x),a(x')} T = 0=[a(x)*,a(x')%, (2.1) 

[a(x),a(x')*] T = 5(x-x'), (2.2) 

where [•,•]_ denotes the commutator, [•,•]+ the anticommutater, and 5(x — x') stands 
for the Kronecker 5 in case of Z 3 . The momentum space operators are distinguished 
only through their argument and given by 

a(k) = (2tt)- 3 / 2 I dxe~ ik - x a{x) , keR 3 , (2.3) 

a(k) = e-^*a{x) , k G [~ i] 3 = T 3 . (2.4) 
(i) Independent particles with random impurities (Anderson model). Particles hop 
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to nearest neighbor sites and are subject to a static potential V(x). Then the 
hamiltonian reads 

H = ~^2 a(x)*(Aa)(x) + V?Y1 V(x)a(x)*a(x) . (2.5) 

x€I? x<=Z 3 

Here A is the lattice Laplacian and {V(x),x G Z 3 } is a collection of independent 
random variables. The random potential mimics the interaction. The kinetic equa- 
tion will be linear, however. 

(ii) Electron-phonon systems. a(x) is the Fermi field of the electrons and c(x) is the 
Bose field of the phonons. The hamiltonian is given by 

H = -- I dxa{x)*Aa{x) + [ dku(k)c(k)*c(k) 

2 JR3 J-^3 

+Ve dx dx'a*(x)a(x)g(x-x')(c(x') + c(x')*) (2.6) 

JR3 Jr3 

with uj the phonon dispersion and g a suitable coupling function. The number of 
electrons is conserved. Assuming the phonon system to be infinitely extended and 
in thermal equilibrium, the dynamics of a single electron is governed by a linear 
kinetic equation. 

(in) Weakly interacting fermions and bosons. If the pair potential is denoted by V, 
the hamiltonian reads 

H = ~f dxa{x)*Aa{x) + / dxdx' a{x)* a(x'YV {x - x')a{x)a(x') . (2.7) 

(iv) Phonons. Phonons are bosons. For effective single band models the free part is 
of the form 

H= I dkco(k)a(k)*a(k) . (2.8) 
./ • ■ 

uj is the dispersion relation, uj 2 is the Fourier transform of a rapidly decaying 
lattice function. For an acoustic band u(k) = \k\ for small k, while for optical 
bands uj(k) > uj > for all k G T 3 . Expanding around the harmonic potentials 
yields cubic and quartic nonlinearities. For example, a quartic on-site potential 
would be A XLgz 3 ^) 4 , Qx the scalar displacement at x away from the equilibrium 
position, while cubic nonlinear springs correspond to | J2 X j/gz 3 a ( x ~ v)(Qx ~ Qy) 3 - 
The respective interaction hamiltonians are 

r 4 
#14 = A / dhdhdk^k^ih + h+h+k^Yli^ik^y^^aik^+ai-kj)*) (2.9) 
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and 



#13 — q / dkidk 2 dk 3 5(ki + k 2 + k 3 ) 

3 J T 9 



<*{x) J! {2u(k j )y 1/2 (e i2 ^ - 1) (a(%) + a(-%)*) . (2.10) 
xez 3 i=i 

Ho+yfeHi 3 is not bounded from below, a defect not yet seen on the kinetic time scale. 
For phonon systems the particle number is not conserved. However, it is possible 
that on the kinetic level only processes with number conservation are allowed. 

3 Initial state 

As explained in Section [TJ for the spatially homogeneous case the initial state is 
translation invariant, quasifree, and invariant under the H dynamics. Therefore 

(o(jfc)) = = (a(k)*) and 

(a(k')*a(k)) = W(k)5(k - k') , (3.1) 

(a(k')*a(k)*) = , (a{k')a{k)) = . (3.2) 

Here W(k) > and for fermions in addition W(k) < 1. All other moments are 
determined through the two-point function, i.e. 

m m 
3=1 3=1 

with all other mononomials having vanishing expectation. Such a state has infinite 
energy and is not supported on Fock space. 

To model the spatially inhomogeneous situation, we give ourselves a spatially 
dependent Wigner function W(r, k), W(r, k) > 0, W(r, k) < 1 for fermions, and of 
rapid decay, in particular J dr J dkW(r, k) < oo. Secondly we choose a sequence of 
trace class operators, R £ , e > 0, on J 7 , R e > and R e < 1 for fermions. Let R e (x, y) 
be the corresponding integral kernel. Then the sequence of initial states, (-} £ , e > 0, 
still satisfies (EO) . (ED, and (E3D with 

(a(xya(y)) £ = R e (x,y). (3.4) 

The Wigner function W appears through the limit e —>■ by requiring 

limiFfe-V + \,e-\- -rj) = [ dke ik ^W{r,k) . (3.5) 

In particular, for the number of particles, N, one has 

(N) £ ~e- 3 . (3.6) 
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(13. 51) implies that at r, on the kinetic scale, in the limit e —>■ the state is 
quasifree with correlator W(r, k)S(k — k') (r is fixed here). Spatial averages on the 
kinetic scale are self-averaging (law of large numbers). To see this, let A be a local 
observable and r x the shift by x. Spatial average on the kinetic scale means 

e 3 dxx\{ex)r x A = n £ (A, A) (3.7) 
</r 3 

with xa the indicator function of a fixed box A. Then 

limn £ (A,A)= / drxA{r)(A) w{r ,) (3.8) 

in the sense that the variance of n £ (A,A) vanishes as e — > 0. (A)w(r,-) is here the 
average of A in the translation invariant, quasifree state with correlator 

(a(tf)*a(k)) w(fY ) = W{r, k)S{k - k') . (3.9) 



4 Kinetic limit 

Let U £ (t) = exp[i(H + y/eHi)t] be the unitary group generated by H. (In brackets 
we note that in all models H is self-adjoint with domain D(Hq) under reasonable 
assumptions on the coefficients, except for H + y/eH^ which has to be stabilized 
by adding a small quartic potential.) As initial state we choose the sequence (-) £ of 
quasifree states from Sectional It is time evolved to (-) e (t) according to (A) e (t) = 
{U £ (t)* All 6 (t)) £ for all A £ -B(jF). To prove the kinetic limit firstly means to 
establish the existence of 

lim^OrV - ^)*a(£- 1 r + -r])) £ (e-H) = [ dke ik ^W{r, k, t) . (4.1) 

e^O 2 2 J R 3 

(For model (i) one averages also over the random potential. For model (ii) one 
considers only a single electron and averages over the thermal state of the phonons.) 
Clearly an issue is the sense in which the limit (14. 1R is intended. Usually one proves 
weak convergence, but pointwise convergence could hold. More ambitiously one can 
study higher moments which should follow the pattern explained in Section [3j Only 
the Wigner function W(r, k) is to be replaced by the time-evolved Wigner function 
W(r, k,t). A further item would be the law of large numbers and its central limit 
corrections. 

Secondly the time evolution of the Wigner function W(r, k, t) must be governed 
by an autonomous transport equation. If Hi = 0, such a property and the limit 
(14. ip are easily established with the result 

%-W(r, k, t) + V k oo{k) ■ V r W(r, k,t)=0. (4.2) 
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If Hi 0, the transport equation (14.21) acquires a collision term as 



-W(r, k, t) + V k u(k) ■ V r W(r, k, t) = C(W) (r, k, t) . (4.3) 

The collision operator C is local in space-time and we only have to display its action 
on the momentum variable. Of course C depends on the particular model. Let us 
give only two examples. For an electron in a random potential one obtains 

C(W)(k) = 2tt / dk'd(k- k')5(uu(k) -u(k'))(W(k') -W(k)). (4.4) 

Here u is the dispersion relation for the electron and d is the Fourier transform 
of the correlator of the random potential, $(x — x') = (V(x)V(x')), (V(x)) = 0. 
The collisions uniformize the momentum over the shell of constant energy. Thus 
with this collision term Eq. (14.31) has the property that as t — > oo the distribution 
function of the electron converges to W(k) = S(u(k) — E), given the initial data 
have kinetic energy E. 

Our second example is the Uehling-Uhlenbeck equation for weakly interacting 
quantum fluids [5], example (Hi), which for fermions was first derived by Nordheim 
[6]. The collision term reads 

C(W)(k 1 ) = J dk 2 dk 3 dk i $(k 1 ,k 2 \k 3 ,h) 

jW 3 w 4 (i + ewi)(i + ow 2 ) - WxW 2 {\ + ew 3 )(i + ew A )] , (4.5) 

Wj = W(kj), with 9 = 1 for bosons and 9 = — 1 for fermions. The collision rate 
encodes energy-momentum conservation and the interaction potential in the Born 
approximation as 

HK k 2 \k 3} h) = \V(ki - k 3 ) + 9V(k 2 - k 3 )\ 2 5(kx + h-h- h) 

5(^1 + ^2-^3-^4), ujj = -k 2 j . (4.6) 

Note that (14. 5 p is in fact cubic. The added term makes symmetry of the collision 
operator more transparent. Phonon collision operators are listed in [7]. They were 
first derived by Peierls [S]. 

One part of the limit theorem must ensure the existence of a unique solution 
to the transport equation ( 14. 3p . For a linear collision operator tools are available. 
For a nonlinear collision operator local in time existence can always be established. 
To have a solution global in time is a difficult issue. Since collisions are at some 
particular spatial location, the nonlinearity is well-defined only if sup r k W{r, k, t) is 
bounded, a property one does not know how to establish in generality. An exception 
are fermions governed by (14.51) with 9 = — 1. Then the solution must satisfy < 
W(r, k, t) < 1 for all t, see [9] for a complete discussion. 
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5 Mathematical results 



For a classical system of hard spheres at low density all facets of the kinetic limit 
are proved under the sole restriction of kinetically short times, i.e. \t\ < to = 1/5-th 
of the mean free time |10j . see also [HI [T2]. The issue under discussion is to carry 
out a corresponding program for weakly nonlinear wave equations, resp. linear wave 
equations with weak disorder. 

The zeroth order step is to have a convincing formal derivation and thereby to 
determine the collision operator. Such formal derivations abound in the physics lit- 
erature. They span from a direct application of Fermi's golden rule to sophisticated 
diagrammatic expansions. To me the most convincing argument is a truncation 
procedure based on the assumption that the state at microscopic time e~H is ap- 
proximately quasifree. The details can be found in [T3], see also [H], for weakly 
interacting quantum fluids and in [7] for weakly interacting phonons. 

Up to now, any proof is based on the time-dependent perturbation theory, which 
means to expand 



Inserting in (U (t)* a(y)* a(x)U (t)) e and working out the average over the initial state 
(resp. over the random potential) this expectation is then represented as a sum of 
high-dimensional oscillatory integrals (= diagrams). Very crudely the next two steps 
consist in dividing between 

• leading diagrams. Their e — > limit does not vanish. 

• subleading diagrams. Their e — > limit vanishes. 
Since (15. ip is an infinite series, one needs on top 

• a uniform bound on the series. 

• the leading diagrams at e = sum to the time-dependent solution of the 
transport equation. 

There are variations in the way how to arrange the perturbation series. The 
expansion of (U(t)*a(y)*a(x)U(t)) e yields standard Feynman diagrams. Particle 
conservation is exploited more efficiently by considering the hierarchy of n-particle 
Wigner functions, n = 1, 2, . . . [T31 [15J, [16] . For a particle in a random potential the 
most powerful method is to use directly (15. ip [T7j . 

A second basic choice is between M 3 and 1? as position space. M. 3 has the 
advantage that for phase factors of the form e~ lk 1 the oscillatory integrals are easily 
estimated. On the other hand the large k behavior needs extra efforts. For Z 3 the 
free propagator is less explicit, but momentum space is T 3 , hence bounded. 




-iHoti 



(5.1) 



8 



For weakly interacting systems the perturbation series diverges on the kinetic 
time scale. This leaves one with analysing individual diagrams. For a translation 
invariant Fermi gas all terms up to order (V{) 2 are studied [I] and the approximation 
W(k, t) = W(k)+tC(W)(k) is proved with C of (14. 5p for 9 = — 1 and integrations over 
T 3 , uj(k) = ^2j =1 (l — cos y). For bosons and fermions moving in M. 3 the same result 
is obtained for spatially inhomogeneous initial data [15]. The next set of results 
concerns the analysis of a certain subclass of diagrams, e.g. all leading diagrams 
[T4"] or the identity permutation in (I3.3P which corresponds to Boltzmann statistics 
[T8] . All these results strongly support the kinetic picture. Spatial dimension has to 
be larger than three, which roughly comes from conditions on oscillatory integrals 
as, e.g., 

POO p 

/ dt\ / dke~ lw{k)t \ < oo. (5.2) 
Jo Ji d 

For the motion in a random potential the situation is more favorable. Assume 
a Gaussian random potential. Then only even moments contribute. For the 2n- 
th term of the perturbation series this yields a bound as ^e n (e~H) n c n ((2n)\/2 n n\), 
where the term in the second round brackets comes from the number of pairings 
for the Gaussian random potential. Thus the perturbation series has a finite radius 
of convergence, to, on the kinetic scale. With this restriction the kinetic program 
can be carried through. The spatially homogeneous case is dealt with in [19], which 
is substantially improved and extended in [20], see also the discussion in [21J. In 
1997 L. Erdos and H.-T. Yau [22] announced a truly remarkable break through. 
By a suitable truncation of the perturbation series, the remainder being controlled 
by unitarity, they succeed to prove the kinetic limit for arbitrary t - With their 
technique other cases can be handled as well, to list: low density scatterers [23] , 
non-Gaussian but approximately independent random potentials [21], tight binding 
dynamics [22], coupling to a phonon bath [26J, and mass disordered lattice dynamics 
[27] . Chen proves self-averaging [2B] and the link to spectral properties of H |25j . 

Staying in the context of the motion of an electron in a random potential, one 
might wonder what happens beyond the kinetic time scale. Mathematically this 
means to consider times of order e~( 1+a H in microscopic units with a > in the 
limit e — > 0. This implies very long kinetic times at which the solution to the 
transport equation is already of the form 

W(r,k,t) = p(r,t)5(E -oo(k)) (5.3) 

with E determined by the initial conditions. The mass density p has a slow variation 
and is governed by the diffusion equation 

^p(r,t) = D Ua (E)A r p(r,t) , (5.4) 

where D^ in (E) is the energy dependent diffusion coefficient of the kinetic equation 

JV(r, k, t) + V k uo{k) ■ V r W(r, k, t) = C(W) (r, k, t) (5.5) 
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with C the collision operator from (14.41) . If a is large, the diffusion coefficient in ( 15.41) 
should be close to the diffusion coefficient as computed for fixed small e, which is 
definitely different from D kin (E). Thus there must be a critical a c up to which (15. 3p 
together with (15.41) is valid. The theoretical prediction is a c = 1. Erdos, Salmhofer, 
and Yau [291 EQ] prove (fO]l . (l5^j) for < a < lCT 3 . 
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